Self-consistent theory of capillary-gravity-wave generation by small moving objects 
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We investigate theoretically the onset of capillary-gravity waves created by a small object moving 
at the water-air interface. It is well established that, for straight uniform motion, no steady waves 
appear at velocities below the minimum phase velocity c m i n = 23cm/s. At higher velocities the 
emission of capillary-gravity waves creates an additional drag force. The behavior of this force 
near the critical velocity is still poorly understood. A linear response theory where the object is 
replaced by an effective pressure source predicts a singular behavior for the wave drag. However, 
experimental data tends to indicate a more continuous transition. In this article, we show that a 
proper treatment of the flow equations around the obstacle can regularize wave emission, even in 
the linear wave approximation, thereby ensuring a continuous behavior of the drag force. 

PACS numbers: 47.35.-i, 68.03.-g 
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I. INTRODUCTION 

An object moving uniformly in an incompressible liq- 
uid experiences a drag force that can have several physi- 
cal origins: viscous drag, hydrodynamic interaction with 
close- by boundaries, or the force due to the emission of 
waves. The waves, on which we focus in the present 
paper, appear when the object moves in the vicinity of 
a deformable surface such as an air-liquid interface 
They carry away momentum from the object which is 
sensed as the wave drag of the moving object. The type 
of waves we expect on an air-liquid interface are cap- 
illary gravity waves 0, @|. Their dispersion relation in 
an unbounded inviscid liquid of infinite depth is well 
known to be ui 2 — gk + jk 3 /p. It relates the oscilla- 
tion frequency u to the wave number k and depends on 
the gravity constant g, fluid density p and on the sur- 
face tension 7. The wave velocity c(k) = co/k is read- 
ily obtained as c(fc) = (g/k + jk/ ' p) x l 2 . The dispersive 
nature of capillary-gravity waves creates a complicated 
wave pattern around a moving object, yielding a finite 
wave drag. In naval design the wave drag is an important 
source of resistance, which stimulated the development of 
approximate theoretical methods d-Q. These methods 
are valid only for objects larger than the capillary length 
— \Jll{P9) @- The case of objects of extension com- 
parable to k _1 has been overlooked for a long time in the 
literature, but has attracted strong interest in the context 
of insect locomotion on water surfaces [IM1- In partic- 
ular, some insect species (for example wiggling beetles) 
may take advantage of the generation of capillary-gravity 
waves for echo-location purposes OH!)]. In particular, 
recent observation of the behavior of Gyrinidae suggest 
that they select their swimming speed by minimizing the 
sum of wave and viscous drags (I6l4l8( j . 

The first theoretical calculation in this regime pre- 
dicted a discontinuity of the wave drag R w at a criti- 
cal velocity given by the minimum of the wave velocity 
Cmin = for capillary gravity waves 0]. For wa- 



ter this evaluates to c m i n ~ 23cm/s. An object moving 
at constant velocity V < c m - ln does not generate steady 
waves, and the wave resistance vanishes. Emission of 
steady waves becomes possible only when V > c m ; n , lead- 
ing to the onset of a finite wave drag. This striking be- 
havior is similar to the well-known Cherenkov radiation 
emitted by charged particles to the onset of wave 
drag for supersonic aircrafts [2(| or to the Zeldovich- 
Starobinsky effect in general relativity The mini- 

mum in the dispersion relation, responsible for this be- 
havior, renders the problem challenging. Two experi- 
ments addressed the problem of the behavior of wave 
resistance at a liquid-air interface. While the disappear- 
ance of wave drag was confirmed for V < c m i n , opposite 
conclusions were reached concerning the existence of a 
discontinuity at the critical velocity V = c m j n . In a first 
experiment by Browaeys et al. [22| the bending of a nar- 
row fiber in contact with the liquid surface was used to 
probe the wave drag. The results evidenced the presence 
of a jump at V = c m - m . However, during these measure- 
ments the contact line between the fluid and the fiber 
was free to move, thus creating an uncontrolled contri- 
bution to the measured force. A second experiment was 
made by Burghelea & Steinberg [23| in which an inge- 
nious feedback system fixed the immersion depth of the 
object. This experiment concluded on a continuous in- 
crease of wave drag around V — c m - m . 

While it was shown recently [24j that the threshold 
V = c m ; n exists only for an object moving at constant 
velocity without acceleration, a theoretical understand- 
ing of the scaling of wave resistance is still missing. The 
theoretical descriptions which are available replace the 
moving object by an external pressure source applied at 
the air-liquid interface. The hydrodynamic problem is 
then reduced to a linear response theory in the pressure 
field which has singularities around V — c m i n if viscosity 
is neglected [9|, [25|, [26J . A self-consistent determination 
of the pressure distribution was attempted in Ref. [261 ]. 
While this theory succeeded in removing the singularity 
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FIG. 1: Geometry of the submerged cylinder and of the free 
surface above. The far-field velocity is indicated by the ar- 
rows V. 



at V = c m ; n , it leads to the somewhat unrealistic predic- 
tion that the applied pressure field vanishes at V = c m ; n . 
An alternative approach by Sun & Keller [13] is based 
on an asymptotic matching technique and works for ve- 
locities much larger than c m i n . In order to develop a 
theory of wave drag valid for V ~ c m i n , we note that the 
linear response theory is successful at reproducing the 
wave pattern created in the experiments even for veloc- 
ities very close to c min . This suggests that it is possible 
to understand the behavior of R w in a theory where lin- 
ear capillary-gravity waves are coupled to an accurate 
hydrodynamic description of the flow around the moving 
object. This theory is developed in the present paper. 

In the situation we have in mind, the perturbation of 
the free surface does not come from an external pressure 
distribution as it was the case in Ref. @. Also, we avoid 
the difficulties arising from an immersed perturbing nee- 
dle of Ref. [26] . Instead, we here use a completely sub- 
merged object, which is, for reasons of a dimensionality 
reduction, a cylinder. It has radius R and is at depth H 
below the free surface, as is depicted in Fig. Q] Notice 
that the dimension has implications on the nature of the 
transition near c m j n which we want to describe and might 
limit comparisons with experimental data. 

The liquid flows only perpendicularly to the axis of the 
cylinder, and we require only x and z coordinates. This 
system has been analyzed by Lamb Q without taking 
into account the mutual interaction between the two per- 
turbations, namely the one created by the moving object 
and the one by the emitted waves. He limited his discus- 
sion to particles larger than the capillary length by 
neglecting the capillary contribution in the dispersion. 
We here regard the other case of small objects, where 
it becomes necessary to treat the mutual interaction be- 
tween the two perturbations correctly. 

Generally speaking, the larger the ratio H/R, the bet- 
ter work the two theories presented below in Sees. ILT1 
andlVl The numerical example in Sec.|lV]for H/R = 5/3, 
however, shows that the theories are valid already for 
cylinders quite close to the surface. This observation is 
the reason why we are convinced that the completely sub- 



merged cylinder presents a useful approach to understand 
the experimental setup in which a partially immersed ob- 
ject was used. 



II. FLOW EQUATIONS 

We describe the velocity field in the reference frame 
of the cylinder. The far-field velocity is thus a non-zero 
uniform velocity V. The velocity field is assumed to be 
irrotational, given by a velocity potential, where we im- 
mediately isolate the far-field velocity and use only the 
potential 4>(x 7 z) of the perturbations due to the obstacle 
and to the waves, 

v(x, z) = V + V(j)(x, z). (1) 

The incompressibility condition for the fluid reads 

= 0. (2) 

This equation is complemented by boundary condi- 
tions at the air-liquid interface and at the surface of the 
cylinder. The kinematic boundary condition reads 



n • (V + V<j 



0, 



(3) 



with the normal vector n oriented as in Fig. Q] For the 
dcformable air-liquid interface the boundary condition is 

pV 2 d 2 x <t> + pgd z (j) - ^d 2 x d z <t> = (at z = 0). (4) 

This boundary condition is obtained as follows. Using 
a height profile C{ x )i the kinematic boundary condition 
becomes for a nearly flat interface (d x ( <C I ), 



- d x {(V + d x <t>) = 0. 



(5) 



Additionally, we have Laplace's law and Bernoulli's equa- 
tion, 



-p(x,C(x))+p ah '- 7 ^ C , 



(6) 



|v 2 (x,C(a;)) + pgC(x) + p(x,((x)) = const, (7) 

with 7 the surface tension and with p a " the pressure 
above the interface, which is assumed to be a constant. 
The curvature of the interface has been linearized. The 
boundary condition (j4|) is obtained by inserting Laplace's 
law into the Bernoulli's equation, linearized in (f>. A 
derivative with respect to x and multiplying with d x (f> 
allows to eliminate the height profile (, using Eq. ([3]). 
Terms which are either quadratic in <p or bilinear in <fi 
and £ are neglected as second-order perturbations. 

The problem we will solve in the following consists of 
the Laplace equation ^ for the unknown <p(x,z), to- 
gether with the two boundary conditions ([3]) and (|4|). 
The flow disturbance must further vanish at very large 
depths: d z <p — > when z — > — oo. These equations can 
not be solved directly by numerical means in real space, 
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FIG. 2: Overview of the variable transformations of the full 
theory in Sees. [TT] and Mil (upper panel), and of the dipolar 
model in Sec. IVl (lower panel). 



such as by a finite element method, since the waves emit- 
ted by the cylinder propagate to infinity whereas any dis- 
cretization is done in a finite domain. Hence, additional 
analytic transformations are needed to cast the problem 
in a more accessible form. An overview of the following 
transformations is provided in Fig. [5J We split the poten- 
tial into two components with different physical origins 
4> = <po + <pi . The term </>o is dominated by the pertur- 
bation stemming from the sphere, and 4>\ is mainly the 
perturbation from the free surface. Of course, both per- 
turbations have to mutually respect the presence of the 
other boundary. In more precise terms, the potential <f>Q 
obeys the Laplace equation with the boundary condition 
n • = —Vc(0) at the surface of the cylinder and 

V^o — > for r — > oo. 

The unknown function Vc{9) = n- (V+V0i) describes 
the flow created by V and the wake at the cylinder sur- 
face (hence Eq. is always verified) . The flow V</>i cre- 
ated by the surface- waves is oscillatory in the x direction 
and decays exponentially with the depth z. These prop- 
erties are naturally captured by the integral representa- 
tion commonly used in theories of surface waves [H, Q , 



i(*.») = ^/ 



dke 1 



>A(k), 



(8) 



where k denotes the wave number. The two unknown po- 
tentials 4>q (x, z) and 4>i(x, z) are now replaced by the two 
unknown functions Vc(0) and A(k). Those have to be de- 
termined self-consistently in order to satisfy the bound- 
ary conditions Eqs. ^ and (g]). This procedure requires 
an expression of (f>o hi terms of Vc , which is given by the 
integral representation 



R 

<f>o(r,9) = / log|re 

7T 



,0 



Re l ° \V c (0')d9' (9) 



which uses the polar coordinates depicted in Fig. Q] No- 
tice that the origin of these polar coordinates coincides 
with the cylinder axis while the origin of the Cartesian 
coordinates x, z is located at the free interface above 
the obstacle. Relation (j9]) is either obtained by expand- 
ing Vciff) as a Fourier series and by solving the Laplace 



equation for each term. An alternative way is to use the 
full Green function of the Neumann problem in the cir- 
cle [IM Hil . The logarithmic integral kernel appears as 
a consequence the fundamental solution of the Laplace 
equation in two dimensions. Notice that the part of the 
Green function which enforces the Neumann boundary 
condition is also logarithmic. 

The boundary condition @ leads to a first relation be- 
tween Vc{0) and A{k). Since it is invariant under trans- 
lations in the x direction, it takes a simple form in Fourier 
space, 

m = -^ u Jl + Zt vM with ( io) 



|fc| uj 2 (k)-V 2 k 2 
V z (k) := I dxe- tkx d z Mx,0) 



(11) 



The expression of V z (k) as a function of Vc(0) is readily 
obtained by inserting Eq. into the above Eq. (fTTl) . 
The integral with respect to x can be evaluated analyti- 
cally, which yields 



v z (k) = -rJ 



e -\k\(H-R sin(e') e -ikR cos(e') ,q/\ ^qI 



The exponential factors arise from the Fourier transform 
of the derivatives of the logarithmic kernel in Eq. © . Wc 
have now expressed the unknown function A(k) in terms 
of the other, Vc(ff). The remaining kinematic bound- 
ary condition ^ at the cylinder surface now serves as a 
closed equation for determining Vc(8), 



Vo(B) = n • (V + V<t>x) 



(13) 



e ikx e \k\z w 2^ + y2 fc 2 



2vr ./ dk \k\ oj 2 {k)-V 2 k 2 



V z (k) 



(14) 



By injecting V z from Eq. dl^j) into the last integral and 
by changing the order of integration, Vc{&) is found to 
satisfy an integral equation of Fredholm's second kind 



V c (ff) = Vcos9+— K(8,6')V c (9')d9', (15) 

Z7T 



with the kernel function 



K(6, Q') = R J [i sgn(fc) cos 9 + sin 9} 



w 2 {k) + V 2 k 2 



L0 2 (k)-V 2 k 2 
^ ^ikR(cos 6 — cos \k\ (2H — _Rsin 6 — Rsin Q') (16) 

In this integral, singular pole contributions appear 
when the denominator vanishes, u! 2 (k) — V 2 k 2 = 0. This 
happens only for V > c m in- The two positive solutions 
of this equation, which we call k g and k c (k c > k g ) cor- 
respond to gravity waves and capillary waves, respec- 
tively. In the far-field regimes, they dominate the flow 
and thus correspond to the inverse wavelengths found be- 
fore and behind the obstacle. Directly above the cylinder 
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no unique wavelength can be identified. In order to en- 
sure that waves only leave the object and are not coming 
back from infinity, we introduce an infinitesimal imagi- 
nary part into the dispersion relation. This term can be 
understood as an infinitesimal viscous term [9j, |26| . With 
the correct choice of its sign, the denominator becomes 
ui 2 (k) — V 2 k 2 — iek. Notice that this choice also ensures 
that capillary waves are emitted to the front and that 
gravity waves rest at the rear of the obstacle. 



III. WAVE DRAG 

Provided that the self-consistent Eq. (|15p has been 
solved to find the function Vc, one can reconstruct the 
flow everywhere in the domain around the obstacle. This 
procedure requires to calculate in turn </>o, A(k), and <f>i, 
using Eqs. (ITU)) , and ©. Once the flow is known, two 
different strategies offer themselves for the calculation of 
the wave drag. Both are indicated in the upper part of 
Fig. [2] The first, which is conceptually simpler, passes 
over the pressure field p{x, z), using Bernoulli's equation 



pv 2 (x, z) 



pgz+p(x,z) 



P V 2 



(17) 



The pressure field can then be integrated around the sur- 
face of the cylinder to find the total drag force. In an in- 
viscid fluid, d'Alembert's theorem ensures that the drag 
caused by the emitted surface waves is the only contri- 
bution to the drag force [H § 11]. The wave drag then 
reads: 



dA{x, z)p(x, z) n. 



(18) 



The second approach to calculate the wave drag makes 
use of the power carried away by the waves. For the 
linear waves used here, the relation between this power 
and the amplitude of the velocity oscillations in the far- 
field regime are well known The oscillations in the 
far- field regime are of the form C,{x) = acos(kx + h). 
The amplitudes a c ,a g , the wave numbers k c ,k g and the 
phases h c , h g are different for capillary and for gravity 
waves. Each energy density of the particular wave de- 
pends quadratically on the amplitude, 



E(k) = ^pc(k) 2 ka 2 . 



(19) 



This energy density may alternatively be expressed in 
terms of the amplitude of the velocities, denoted by v ZiOC , 
which is related to the amplitude a of the deformation 
by f 2 ,oo = Vka. This relation is a consequence of the 
linearized kinematic boundary condition ([5]). It allows 
to determine the amplitude a from the flow, once the 
function Vc{9) is known. The two energy densities allow 
to calculate the power carried away by the waves, from 
which the drag results as Q, 



V-uj'(k c ) 
V 



E(k c ) 



V 



4^E(k g ). 



(20) 
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FIG. 3: (Color online) Numerical solution of Eq. lfl5)l. The 
dashed lines represent the discretisation with linear finite el- 
ements. Parameters: Hk = 0.5, Rk — 0.3, and V = 1.3c m j n . 



The energy of each wave is transported at the group ve- 
locity uj' . 

The two expressions (fT8|) and (|20|) for the wave drag 
are equivalent. Their comparison allows to verify the 
accuracy of the numerical solutions below. 



IV. NUMERICAL RESULTS 

The self consistent equation Eq. (|15p can be solved by 
numerical methods. Our first step is to create a numeri- 
cal table of the values of the kernel K(9, 9'), the integral 
in the definition of the kernel is computed using stan- 
dard numerical routines from the GSL library [30|. It 
provides special discretizations for calculating principal 
values which appear in Eq. (|16[) due to the singular in- 
tegrand. The equation Eq. (1151) is then discretized using 
finite elements. The function Vc{9) is approximated us- 
ing a finite set of N basis functions 4> n (6) . 



N 



(21) 



For simplicity we have chosen a basis of N piecewise lin- 
ear hat functions with periodic boundary conditions on 
the (— 7r, 7r) interval. Fig. [3] shows an example of such an 
approximation for a small number of elements N = 20. 
This expression for Vc(0) is then inserted in Eq. (fl5|) 
and Galerkin's method is used to convert it into a matrix 
equation of size N x N. Namely we multiply Eq. (IT5l) by 
one of the function <j) m {6) and integrate over the angle 9. 
This leads to a linear equation on the coefficients c n : 



E 



i 

2^ 



(cj) m ,K<p n )]c n = V(<p m ,cosd) (22) 



where we have introduced the notation K<j> n (ff) = 
J K(6,8')4> n (0')d9' and defined scalar product between 
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FIG. 4: (Color online) Visualization of the flow velocity 
around the cylinder. The velocities are represented by ar- 
rows the size of which is proportional to the velocity magni- 
tude. The color code is yellow/gray for large magnitudes and 
blue/black of small ones. Parameters as in Fig. 



two arbitrary functions f{9) and g{9) as 

(/,<?):= / f(0)g(0)d8 (23) 



In a last step the system Eq. (|22l) is solved using a LU 
decomposition [3l[. The integrals involved in the cal- 
culation of the scalar products are determined numeri- 
cally [3(| • It is possible to check the convergence of our 
numerical procedure by inserting the obtained approxi- 
mation for Vc(0) into Eq. (fl"5|) . The degree of accuracy 
can then be estimated from the difference between the 
left and right hand side of Eq. (|T5|) . For the typical num- 
ber of elements we use in our simulations N ~ 110 the 
relative difference is of the order of 10~ 3 . 

Figure [3] presents the solution of Eq. (|15[) for the ge- 
ometrical parameters Hn = 0.5, Rk = 0.3 and a flow 
velocity of V = 1.3c mm . It is compared with the source 
term of the integral equation V cos 9 which results if the 
contribution from the surface waves is neglected. As can 
be seen, the function Vc{6) is strongly modified by the 
self-consistent interaction between the cylinder and the 
emitted capillary-gravity waves. 

Once the function Vc{9) has been determined we can 
calculate the flow in all the space z < around the cylin- 
der following the steps described in section II. The veloc- 
ity field around the cylinder for the parameters of Fig. [3] 
is depicted in Fig. 2J Notice that the flow obeys the ki- 
netic boundary condition on the sphere given by Eq. (|3]) 
which confirms our numerical procedure. 

Using the results from Sec. IIII1 we can now calculate 
the wave drag as a function of the externally applied flow 
velocity V. Figure [5] shows this function, as determined 
numerically from both, the direct integration of the pres- 
sure field around the sphere (|18[) . and from the energy 
balance in Eq. (|2"0)l . Their difference in Fig. [5] is smaller 



1.5 



.5 1 



==0.5 - 



Linear response 
Self consistent model 
Dipolar model 




1.5 



V/ Cmin 



2.5 



FIG. 5: (Color online) Wave drag as a function of the far-field 
velocity V for the parameters for Hn — 0.5 and Rk — 0.3. 
The results of three different models are presented: the dashed 
curve corresponds to a linear response calculation where Vc = 
VcosO, the red (gray) curve corresponds to a self-consistent 
calculation using Eq. {15} and the black curve corresponds to 
the dipolar model explained in section IV. 



than the line width. At velocities V close to the critical 
value c m in, the drag vanishes continuously. This behav- 
ior is qualitatively different from that found from linear 
response theory with an external pressure field as per- 
turbation, as it predicts a divergence at c m ; n @. Our 
approach reduces to this linear response theory if, in- 
stead of solving the self consistent Eq. (fT5l , Vc (9) is set 
to the result V cos 9 for a cylinder moving in a quiescent 
inviscid liquid without surface. The wave drag obtained 
in this case diverges close to V — c m i n as shown by the 
dashed curve on Fig. [5] At large velocities, the wave drag 
determined from the self consistent model curve rejoins 
the curve from linear response. 

The observation in Fig. [51 that the drag is a continuous 
function at the critical velocity, constitutes the central 
result of the present work. The regularization is achieved 
by the self-consistent treatment of the emitted waves, 
which takes into account the mutual interaction between 
the perturbations created by the cylinder and the waves, 
respectively. The physical origin of the regularization 
can be understood with a simpler model which does not 
require the full solution of the self-consistent equation. 
This model will be detailed in the following section, where 
we derive a square-root scaling for the wave drag close 

to C m j n . 



V. DIPOLAR MODEL 

We now treat the interaction between the wake and the 
cylinder in an approximate manner. Instead of enforcing 
the exact kinematic boundary condition at all points 
on the cylinder, we impose it only as an averaged con- 
straint. This approach assumes that the flow vi = V0i, 
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created by the waves, is homogeneous around the cylin- 
der. The response of the cylinder is now reduced to a 
simple dipolar potential which describes the response to 
a yet unknown uniform flow (V +Vi), 



i? 2 (V + vi)- Vlogr. 



(24) 



The name "dipolar" arises from the analogy with the 
potential of a dipole in two-dimensional electrostatics. 
This approximation reduces the complexity of the equa- 
tion (J9j> . Figure [2] summarizes the necessary transfor- 
mations in the same way as it was done for the sections 
above. Instead of determining a whole function 0o we 
must now only find two parameters vi yX and v\ tZ . With 
these parameters, the function Vc{0) from Sec. [TT] as- 
sumes the form Vc (0) — (V+Vi iX ) cos 0+v\^ z sm.0, which 
can be used immediately to calculate the unknown func- 
tion A{k). The self-consistent equation (TT5|) for Vc is 
thus replaced by a linear set of two equation for the pa- 
rameters. The solution is 



V + v hx \ = V_ 

with the shortcuts 



I 2 



h = 
h = 



dk 

— ttR e 
2tt 

dk 

2tt' 



2-2H \k\ 



7rR 2 e~ 2H ^ik 



LJ 2 {k) + V 2 k 2 

w 2 (k) 

UJ 2 {k) - 



- V 2 k 2 ' 
V 2 k 2 



Lo 2 (k)-V 2 k 2 ' 



(25) 

(26) 
(27) 



In order to obtain the wave drag R w from this result, 
we follow the same steps as in Sec. IIII1 This time, we 
can do them analytically, not only numerically, thanks 
to the simple form of the potential <f>Q. In this treatment, 
we prefer the energy balance argument instead of the 
pressure integral. In the far- field, x — > ±oo, we find the 
wave amplitude as follows, 



= 2-kR' 



kVe~ H ^ 



1 



k(fc) - v\ ^(i - hy + q 



c(k) = V (28) 



In the course of this calculation, we pass from (f>o to 
A(k), using Eqs. (ITU1) and (TiT]) . In the far field regime 
the only contribution to the flow arises from the waves. 
The flow potential </>i created by the waves is given by a 
Fourier-Laplace transform of A(k) (see Eq. (jHJ). In the 
asymptotic regime \x\ — ¥ oo we can keep only the con- 
tributions from the emitted waves of fixed wave numbers 
k = k c ,k g . They appear from the delta function contri- 
bution of the integration around the poles of A(k), using 
the identity [H p. 481]: 

1 



l^o uj 2 (k) - V 2 k 2 - iek 

P.v. — 1 + i7rsgn(fc) S(uj 2 {k) - V 2 k 2 ), 



u 2 {k) - V 2 k 2 



(29) 



where P.v. denotes the Cauchy principal value and the 
delta function can be transformed as: 



5{u 2 {k)-V 2 k 2 ) = £ 



ki:c(ki)=V 



2V\h\ \u'(ki)-V\ 



(30) 



Now, that we have an expression for the deformation 
amplitude a, we can find the wave drag using Eqs. (|19[) 
and ([2U)l . The auxiliary integrals I\ and I2 are calculated 
numerically. The result are displayed in Fig. [5] and show 
that the dipolar model reproduces the behavior of the 
exact solution from Sec. [TTJ The accuracy of the dipo- 
lar model can be understood because the radius of the 
cylinder is small compared to the capillary wave length 
R = O.'Sk^ 1 . Hence the flow is reasonably uniform on this 
scale. A good agreement is found when V —> c m - ln because 
the emitted wavelengths tend both to the same value 
which is larger than the object size. At higher velocities, 
where the emitted wavelengths differ, the agreement de- 
creases because the shorter one may become smaller than 
the object size. 

Our aim now is to find the square-root behavior of the 
drag force, R w cx \JV — c m i n in the frame of the dipolar 
model. In order to do so, we analyze the scaling of the 
wave amplitude in Eq. (|28|) in the vicinity of c m - ln . We 
then use Eqs. p^|) and (j2"0)) to establish the scaling for 
the wave drag. All scaling behavior near c m i n stems from 
the scaling of \u>'(k) — V\ where k obeys c(k) = V. Near 
k = k the phase velocity shows a quadratic minimum and 
the selected wave numbers are: 



V 



(31) 



The difference \oj'{k) — V\ = |w'(fe) — c(k)\ is thus deter- 
mined by the linear contribution from the group velocity 
which is given by: 



(32) 



Note that the result is the same for both capillary and 
gravity waves. 

The scaling of \ui'(k)— V\ is taken over to the amplitude 
a, to I2 and to the wave drag R w . The term I\ becomes 
constant in the limit V — > c m i n , while the term I2 di- 
verges. Indeed 1\ involves only the principal value term 
from Eq. (|2"5)l while I2 depends on the delta function part 
that scales as l/\uj'(k) — V\. After all these preparations, 
we are now ready to provide the scaling of the wave drag: 



R,, 



P C n 



V 



(33) 



This expression is derived assuming R <^ H k _1 , and 
it holds for velocities very close to c m j n , where I2 3> 1 
(this occurs as soon as V — c m i n <C (n 2 R 2 ) 2 ). 

While our approach is exact for a linear capillary- 
gravity waves non linear corrections become relevant 
when the flat interface condition d x C, <C 1 is not verified 
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anymore [33[. We have checked numerically the valid- 
ity of this assumption for the parameters of Fig. [5] It 
is clearly verified at large velocities but as V approaches 
c m i„ the typical d x (^ increases and reaches a regime where 
d x C, — 1 for V ~ c m i n . In this range, where we find 
moderate nonlinearities of the order one, we expect that 
they will change only the quantitative predictions of the 
self-consistent theory, without introducing qualitatively 
new features. This is consistent with the observation 
that nonlinear waves (for e.g. solitons) are not seen in 
the experiments close to c m in and that form of the gen- 
erated wave pattern is well described in linearized the- 
ory [23|]. Hence although nonlinearities should be taken 
into account for a complete hydrodynamic theory, the 
regularization of the wave drag near V ~ c m ; n is already 
present in the linear self-consistent theory. Notice that 
approaches that include only nonlinear effects without a 
self-consistent treatment fail to produce stable solutions 
near V = c min [34j . 

Recently, Moisy, Rabaud, and Salsac [35[ developed an 
experimental technique that allows to recover the height 
profile created by a moving disturbance. This technique 
is a promising candidate for the verification of our theo- 
retical predictions on the onset of the wave drag and the 
implied wave pattern close to c min . 

VI. CONCLUSION 

We have addressed the behavior of the wave drag for 
objects of extension smaller and close to the capillary 
length moving at a speed close to the critical velocity 
Cmin, which is given by the minimum of the wave veloc- 
ity for capillary gravity waves. It is known that theo- 
ries where the object is modeled by an external pressure 
source lead to singular wave drag behavior at V — c mra - 
In this article, we show that even for linear capillary- 
gravity waves, this singularity can be removed, even in 
the approximation of an inviscid fluid, by treating the 
boundary conditions at the object interface in an exact 
way. For this purpose, we treat the wave emission prob- 
lem by a cylinder submerged near the liquid-air interface 
under an external flow. We derive a self-consistent inte- 
gral equation describing the flow velocity at the cylinder 
interface. This equation is solved numerically with a fi- 



nite elements method, which allows us to reconstruct the 
flow in the entire space around the cylinder and to deter- 
mine the wave drag on the cylinder. In addition to the 
numerical solution, we propose a simple approximation 
valid for cylinder diameters smaller than the wavelength, 
where the interaction between the waves and the cylin- 
der is treated in a dipolar approximation. In this case 
it is possible to make analytic estimates showing that 

Our findings explain why a smooth onset of the wave 
drag is observed even if the shape of the wave pattern is 
well described by linear response. The validity of our the- 
ory is limited by the validity of the flat interface approx- 
imation. In our simulations, this conditions is not well 
verified for velocities around c mm . In principle, nonlinear 
corrections should therefore become relevant. However, 
we think that the qualitative behavior is already cap- 
tured by our self-consistent theory, even if the inclusion 
of nonlinear effects would be required for quantitative 
predictions. 

One might be tempted to compare the scaling R w cx 
y/V — Cmin, which we found in the dipolar model, with 
the square-root fit done by Burghelea & Steinberg [23[ 
in their Fig. 18. This apparent agreement between the- 
ory and experiment must be taken with care. Our results 
are valid for a long cylindrical obstacle, while the exper- 
iment was done with a spherical object. Moreover, the 
experimental data does not allow to determine unam- 
biguously a square-root scaling at the onset of the wave 
drag. Hence, more experiments and a three-dimensional 
theory will be required to establish the exact scaling at 
the transition. The main point however, namely the fact 
that our theory is able to recover the continuous onset 
of the wave drag at Cmin, coincides with the mentioned 
experiment. Our finding of a continuous drag force might 
serve as an element of understanding the motion of small 
insects on or near water surfaces, since such animals have 
to find a delicate balance between viscous and wave drag. 
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